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ABSTRACT 

We investigate the epicyclic motion of stars escaping from star clusters. Using streak- 
lines, we visualise the path of escaping stars and show how epicyclic motion leads to 
over- and underdensities in tidal tails of star clusters moving on circular and eccentric 
orbits about a galaxy. Additionally, we investigate the effect of the cluster mass on 
the tidal tails, by showing that their structure is better matched when the perturbing 
effect of the cluster mass is included. By adjusting streaklines to results of A^-body 
computations we can accurately and quickly reproduce all observed substructure, es- 
pecially the streaky features often found in simulations which may be interpreted in 
observations as multiple tidal tails. Hence, we can rule out tidal shocks as the origin 
of such substructures. Finally, from the adjusted streakline parameters we can verify 
that for the star clusters we studied escape mainly happens from the tidal radius of 
the cluster, given by xl — {GM/{^'^ — We find, however, that there is 

another limiting radius, the "edge" radius, which gives the smallest radius from which 
a star can escape during one cluster orbit about the galaxy. For eccentric cluster orbits 
the edge radius shrinks with increasing orbital eccentricity (for fixed apocentric dis- 
tance) but is always significantly larger than the respective perigalactic tidal radius. 
In fact, the edge radii of the clusters we investigated, which are extended and tidally 
filling, agree well with their (fitted) King radii, which may indicate a fundamental 
connection between these two quantities. 
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1 INTRODUCTION 

Tidal tails of galactic satellites like dwarf spheroidal 
galaxies and globular clusters are interesting structures 
as they tell the story of the dynamical past of the corre- 
sponding object. The most prominent feature of tidal tails 
is their large extent across their ho s t galaxy (|Yannv et aL 
2003; Grillmair fc JohnsonI l2006l : iGrillmair fc Dionatos 
2006ai ibf). From the shape of extended tidal tails it is 
possible to infer t he orbit o f their progenitor satellite 
jLin fc L vndcn-B elj Il977l: 
19951: iPcfiarrubia et al.l |2005| : 
Fellhaue~et al.l |2006| ). The constraints on the 
orbit can be further improved when radial 



BinnevI l2008l: lEvre fc B inncv' '2009allbl: lOdenkhchen et all 



20091 : IKoposov. Rix fc Hogg„201Q ')" 



Lvnden-Bell fc Lvnden-Belj 
Belokurov et all l2006l : 
satellite 
velocity 

or proper motion da t a of tidal tail stars is available 
|johnston et all 1 19991 : iLaw. Johnston fc Maiewskil l2005l : 
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Tidal tails, however, also contain additional informa- 
tion. Variations of the stellar density along the tails have 
been observed for extended tidal streams such as the Milk y 
Way globular cluster Palomar 5 l|Odenkirchen et al.ll2003l V 
For other Milky Way satellites, such as NGC 288, Will- 
man 1 and NGC 2298, multiple tida l tails have been ob- 
served jLeon. Mevlan fc Combesll200ol : IWillman et"al]|2006l : 
iBalbinot et al.l bOUl ) This substructure should also be 
somehow linked to the dynamical past of the satellite. 

In the common picture, substructure is created by tidal 
variations, such as tidal shocks, which temporarily increas e 
the mass loss rate of the satellite jGnedin fc Ostrikeijll997j V 
Tidal tails which are produced in such a violent way should 
be dynamically hotter and thus broader than tails origi- 
nating from dynamical evaporation. Extended satellites like 
dwarf spheroida l galaxies often show such tidal streams (e.g. 
iGrillmaiil |2006| ) but tidal tails of st ar clusters appear to 
be typically more compact and colder (|Grillmair fc JohnsonI 
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l2006l :[ Grillmair fc Dionatodl2006bl ) suggesting that they are 
mainly being produced by dynamical evaporation. 

In numerical investigations it has been shown that for 
star clusters, tidal variations have to be, in fact, very ex- 
treme in order to cause substructure in their tidal tails 
l|Dehnen et al.l |2004| : iKiipper et all l2010al '). This is due to 
the compactness of star clusters which makes them less vul- 
nerable to tidal stripping. Observed substructure in tidal 
tails of star clusters is, therefore, most likely not due to 
tidal variations. 

On the contrary, it has been shown that substructure 
can even arise in t idal tails of star clusters which move in a 
steady tidal field llCapuzzo Dolcetta, Di Matteo fc Miocchil 



20051: iKiipper MacLeod fc Heggid l2008al : lJust et al.l l2009l : 



Kiipper et al.l l2010al ). This substructure is caused by the 



epicyclic motion performed by a continuous stream of stars 
evaporating from the cluster. It has been shown analytically 
and numerically that this substructure can be related to 
the internal properties of the star cluster (its mass) and its 
external properties (its orbit about the galaxy). Moreover, 
star clusters on eccentric orbits were also found to produce 
tidal ta il substructure in this way, rather than through tidal 
shocks (|Kiipper et allbOlOal ). 

However, clumps in tidal tails may also be produced 
when parts of a stellar stream pass a spiral arm, a gi- 
ant molecular cloud or a massive dark-matter subhalo 
(e.g. Ilbata erai]|2002l : ICarlberg,. 2009, 1 . Therefore, the under- 
standing of the formation of substructure in tidal tails is of 
vital importance when interpreting observations of clumpy 
tidal tails. 

In the present investigation we show with the help of 
streaklines exactly how the epicyclic motion of escaping stars 
leads to substructure in tidal tails of star clusters on circular 
and eccentric orbits (Sec. [2| . Then we compare results from 
A'^-body computations to theoretical streaklines to visualise 
the path of escaping stars and to demonstrate the accuracy 
of the much more rapid method of streaklines (Sec. [3]). In 
the last Section we give a brief summary and discuss our 
results. 



2 THEORY 

Due to energy equipartition, stars in a star cluster perma- 
nently exchange energy. This process is called two-body re- 
laxation and leads to the continuous escape of stars from the 
cluster, since again and again some stars gain energies above 
the escape energy (see e.g. Hoggic & Hut 2003). However, 
in most cases, even if they are energetically unbound, stars 
will stay within the cluster for several dynamical timescales. 
These "potential escapers", are preferentially located at 
large cluster radii where the dynamical timescale is of the 
order of the orbital timescale of the s t ar cluster about the 
Gala,xv |Gnedin. Lee fc Qstriker|[l999l : iFukushige fc Heggid 
I2OOOI : IKiipper et aLlbOlObf T In fact, a flattening of the veloc- 
ity dispersion profile has been observed in the outer parts oi 
Milky Way globula r cluster s like 47 Tucanae (Drukier et al.' 



1998 : IScarpa et aU2007: L ane et al.|[2QQa; Scarpa fc Falom o 
2OI0I : iLane et al.rr2010a| [bl: lscarpa et ahlfeoili ). which mav 



well be interpreted as the contribution from potential esca- 
pers. 

External energy which is added to a cluster, for in- 



stance through tidal shocks, is redistributed within the clus- 
ter in the same way. The additional energy causes more 
stars to escape from the cluster in a given time interval but 
does not immediately influence the structure of t he cluster 
iBaumgardt fc Makindlioosi : IKiipper et al.ll2010bl ). For this 
reason tidal shocks do not inevitably lead to substructure 
in tidal tails as the escape of unbound stars is d elayed with 
respect to the tidal shock jKiipper et a l."2010a''). Note that 
this may be significantly different for extended objects with 
large relaxation times such as dwarf spheroidal galaxies, es- 
pecially when t hey are on very eccentri c orbits about their 
host galaxy fc.f. iPefiarrubia et al.|[2009l ). 

Inside the cluster, each star moves within the gravita- 
tional field of the other cluster stars plus the external tidal 
field of the galaxy. In the cluster reference frame, in which 
the cluster is at rest, the motion of stars is complicated by 
the Coriolis, centrifugal and Euler pseudo forces resulting 
from the cluster's orbita l moti on about the galactic centre 
(see e.g. IChandrasekhaij 1 19421 ). This effective potential of 
the cluster is attractive towards the cluster centre within 
the tidal sphere with radius xl- This radius gives the dis- 
tance to the Lagrange points from the cluster centre along 
the radial direction connecting the galactic centre with the 
cluster centre, i.e. the radius at which the cluster attrac- 
tion equals the effective force of the external field. It is often 
denoted as the tidal radius and can be approximated for 
circular cluster orbits in a Milky Way-like potential as 



XL 



(GM_\ 
\2Q?) 



1/3 



(1) 



where G is the gravitational constant, M is the cluster mass, 
and is the cl uster's angular velocity abo u t the galactic 
centre (see e.g. iKin 3 119621 : lJust et al.|[2009l : IKiipper et al.l 
l2010al V For eccentric cluster orbits the tidal radius can be 
calculated using 



XL = 



GM 



1/3 



(2) 



where $ is the galactic pote ntial and R is the cluster's galac- 
tocentric distance (e.g. Kiip per et al.l[2010al . Due to the de- 
pendence on n, the tidal radius can change significantly be- 
tween peri- and apogalacticon for clusters on eccentric or- 
bits. When a star evaporates from a cluster it will most likely 
escape through one of the two Lagrange points L\ (closer to 
the galactic centre) or L2 (further away from the galactic 
centre), since escape is easiest through these points, finally 
ending up in the leading or trailing tidal tail, respectively. 

Note that the estimate of the limiting radius (Eq. (2} 
was derived for a star instantaneo usly situate d between the 
cluster and the galactic centre (see lKin3ll962l ). Considerable 
subsequent research has shown that stars well outside this 
estimate of the limiting radius may well remain attached 
to the cluster, for e xample those on appropriate retrograde 
orbits l|Read et al ] |2006l ). Hence, we should keep in mind 
that at any time there will be stars outside the radius xl 
which are actually still bound to the cluster, and that there 
are stars escaping from radii smaller than xl- 

Outside the tidal radius, the effective potential is repul- 
sive in some directions such that stars feel a drag away from 
the cluster once they escape from its tidal sphere. Hence, 
the gravitational field of the cluster quickly becomes unim- 
portant for the motion of the escapers. Within the tails the 
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stars will perform an epicyclic motion due to the pseudo 
forces mentioned above (plus the tidal force), in which the 
stars will be periodically accelerated and decelerated while 
moving along the tidal tails. How this epicyclic motion looks 
depends on the escape conditions of the stars. In the case 
where the escape conditions of many escaping stars are sim- 
ilar, their epicyclic motions will lead to an over- and under- 
density pattern. This is a statistical effect, which is due to 
the fact that there is a higher probability for the presence of 
stars within the epicyclic cusps than in between two cusps. 

In the following we briefly review how the epicyclic mo- 
tions of escaping stars depend on the escape conditions. For 
three different star cluster models we then show how the 
orbit influences the epicyclic trajectories of escaping stars. 
Moreover, we demonstrate the influence of the cluster mass 
on these trajectories. Since the stellar trajectories cannot 
be calculated analytically for clusters on eccentric orbits or 
when the cluster mass is taken into account, we use streak- 
lines as often applied in fluid dynamics to visualise the tracks 
of escaping stars (see Sec. 12. 3[) . 



2.1 Epicyclic motion 

In lKiipper. MacLeod fc Heggid (|2008al ) it was shown for star 
clusters on circular orbits that, in a coordinate system in 
which the cluster is at rest (x-axis points towards the galac- 
tic anticentre, y-axis points into the direction of cluster mo- 
tion), the length of an escaper's epicycle along the corre- 
sponding tidal tail (leading or trailing), yc, can be analyti- 
cally related to the tidal radius, xl, oi the cluster by 



yc = 



1 - 



XL, 



(3) 



where k is the epicyclic frequency of the specific galactic 
potential. For clusters in M ilky Way-like tidal fields k ~ 
lAn holds l|Just et al.ll2009l '). such that Eq. [3] reduces to 



yc ^ Stvxl- 



(4) 



In this calculation it is assumed that the escaping star evap- 
orates from one of the two Lagrange points and has exactly 
the same angular velocity, Q, as the cluster at the moment 
of escape. Moreover, the influence of the gravitational at- 
trac tion of the cluste r is neglected in this ansatz. 

lJust et al. I l|2009l ) found a more general solution for the 
length of the epicycles. 



yc 



4tv 13^ 



13 13^ 



(ax + ^) 



(5) 



where /3 = k/Q, Ax is some distance from the cluster centre 
to the point of escape along the radial direction connecting 
cluster centre and galactic centre, and vt is some additional 
velocity along the cluster orbit in the direction of the corre- 
sponding tail. For a Milky Way-like potential this simplifies 
to 



yc 



(6) 



which gives Eq. [4]for Aa; = xl and vt ~ 0. Therefore, if a 
star escapes from a radius larger or smaller than the tidal 
radius, the epicycle gets larger or smaller, respectively. Fur- 
thermore, an additional tangential velocity in the direction 
of the corresponding tail will increase the epicycle length of 



a star, whereas an additional tangential velocity in the op- 
posite direction to the corresponding tail will decrease the 
epicycle length. 

The radial offset and the radial amplitude of the 
epicyclic motion (the x-com ponent of the epi cyclic motion) 
was derived analytically bv I Just et al.l l|2009fl . This part of 
the epicyclic movement depends on the initial radial offset. 
Ax, and on the additional tangential velocity along the tidal 
tails, Vt- It is given by 



(7) 



where we corrected a typo in the original Just et al. equation. 
This amplitude reduces in a Milky Way-like potential to 



2A.+ g± 



Ax -I- 



(8) 



that is if Vt is zero then the escapers oscillate around a radius 
of 2A3; with an ampl i tude of 2 Ax. 

In iKiipper et al.l Hoioi) it was shown with a compre- 
hensive set of A'^-body models that the epicycles of escaping 
stars from clusters on circular orbits are on average a bit 
larger than predicted when assuming Aa; = xl and vt ~ 0. 
In the framework of this theory, this increase in length can 
be due to both Ax > xl or vt > 0. In the following we 
will demonstrate that this discrepancy is, in fact, due to 
the neglected cluster mass. Since there is no simple analyt- 
ical solution to the problem when the cluster mass is taken 
into account, we can only prove this assumption numerically, 
which will be done using streaklines. 

Fu rthermore, it was demonstrated in iKiipper et al.l 
(|2010ah that epicyclic motion of escaping stars also leads to 
over- and underdensities in tidal tails of star clusters mov- 
ing on eccentric orbits. Even though the shape and density 
distribution along the tails gets significantly more compli- 
cated for such clusters, it was shown that the distance of 
the epicyclic overdensities could be related to a mean es- 
cape radius and a mean escape velocity. Due to the periodic 
acceleration and deceleration of the cluster and the tail stars 
within the orbit about the galactic centre, the distance from 
the cluster to the overdensities periodically increases or de- 
creases, respectively. For this problem there is no simple 
analytical solution either, therefore, in the following we are 
also going to visualise the trajectories of escaping stars for 
such cluster orbits using streaklines. 

2.2 Models 

We investigate the epicyclic motion of tail stars for three 
different types of cluster orbits. First of all, we look at the 
circular orbit case, then we look at two eccentric orbits with 
eccentricities, e, of 0.25 and 0.5. The eccentricity is defined 

as 



Ra 



Rn 



Rapo ~h Rpe 



(9) 



where Rapo is the cluster's apogalactic distance and Rperi 
its perigalactic distance. A circular orbit therefore has an 
eccentricity of 0. 

For the integration of t he orbits we use th e Milky Way- 
like potential suggested bv lAUen fc SantillanI (|l99lh . This 
galactic potential, "I>, consists of a central point-mass poten- 
tial given by 
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Figure 1. Orbits of the three A'^-body models which were used 
in this analysis for comparison with the theoretical streaklines. 
One model is on a circular orbit, one has an orbital eccentricity 
of 0.25 and the third model follows an orbit with an eccentricity 
of 0.5. The crosses mark the orbital phases where snapshots were 
taken for analysis (see also Tab.[T]l. 



Table 1. Overview of the snapshots from the three cluster mod- 
els which are investigated in detail. The columns give the orbital 
eccentricity, e (Eq. [9]|, the apogalactic distance, Rapo, the peri- 
galactic distance, Rperi, and the galactocentric distance, Rgc ^ 
of the cluster. The orbital phase, Porti is defined by Eq. 1141 The 
two last columns give the orbital velocity, i>or6i ^'^'^ the angular 
velocity, f2, of the cluster at the time of the snapshot. 
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a lMivamoto fc Nagail l|l975h disk potential given by 



(10) 



(11) 



and a halo potential of the form 
M(R) Ms 



R 1.02a3 



(12) 



1.02 



1 + Was) 



— + ln{l + {R/a,y-''') 



where 



M(R) = 



1 + [R/aa) 



(13) 



The numerical constants Mi, bi, M2, 02, 62, M3 and 03 
are chosen such that the combined potential of the three 
components has a circular velocity of 220 kms^^ at 8.5 kpc 



com ponents has a circular y e 
(see Uuen fc Santillanlll99ll ). 



The apogalactic distance of all three clusters is 8.5 kpc 
(which is a somewhat arbitrary choice). Therefore, the circu- 
lar orbit has a constant orbital velocity, Vorb, of 220 kms~^, 
whereas the eccentric orbits have an orbital velocity of 165 
kms~^ and 110 kms~^ at apogalacticon for an eccentric- 
ity of 0.25 and 0.5, respectively. Their perigalactic distances 
are 5.1 kpc and 2.8 kpc, respectively, and their velocities 
at these orbital phases are 275 kms~^ and 334 kms~^. The 
orbital timescale is 240 Myr for the circular case, 133 Myr 
for e — 0.25 and 113 Myr for e = 0.5. An overview of these 
quantities is given in Tab.[T] In Fig. [T] the integrated orbits 
are shown for 400 Myr. 

The shape and structure of tidal tails for star clusters on 
eccentric or bits is strongly influenced by the orbital phase of 
the cluster l|Kiipper et al ] |2010al : lKiipper7 Mieske fc Kroupal 
L2011a). To illustrate the changes between the different or- 
bital phases, we pick four snapshots from the integrations of 
each of the two eccentric orbits. The snapshots are taken 
in steps of 30 Myr starting from apogalacticon. In the 
next snapshots the clusters are in an orbital phase between 
apogalacticon and perigalacticon in which the clusters are 
accelerated towards the galactic centre. The third snapshots 
are close to perigalacticon, and the last ones are between 
perigalacticon and apogalacticon in which the clusters and 
their tails are decelerated. The points at which those snap- 
shots were taken are illustrated in Fig. [T] as crosses. 

One way to quantify the orbital phase, Porb, of a 
star cluster in its orbit about the galaxy is given in 
iKiipper. Mieske fc Kroupal l|2011al ') by 



porb 



Rgc Rgc — Rper 

I -Rgc I ^a-po ^ Rper 



(14) 



where Rgc is the cluster's current galactocentric radius, and 
Rgc is the time derivative of this radius. In this way, the 
orbital phase is zero at perigalacticon and unity at apogalac- 
ticon (in fact ±1). Moreover, it is negative when the cluster 
is moving to perigalacticon and positive if its moving to 
apogalacticon. The orbital phases of the cluster snapshots 
which we investigate in detail are given in Tab. [T] 



2.3 Streaklines 



In iKtipper. MacLeod fc Heggij l|2008al ') and lJust et all 
(l2009l ) it was shown that the orbits of stars escaping from 
a cluster can, under certain assumptions, be integrated an- 
alytically. In this simplified case a star is released from the 
cluster centre with a certain positive or negative offset along 
the galactocentric radius, and with some negative or positive 
offset velocity with respect to the cluster. One of the main 
assumptions here is that the cluster moves in a circular or- 
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Figure 2. Snapshot of the cluster model on a circular orbit (black dashed lines). Left panel: to produce the three sets of streaklines, 
test particles were released during the last few hundred Myr of this numerical integration from two radial offsets, ±Ax (which we chose 
to be 60 pc in this example), with respect to the cluster orbit, and with three different tangential velocities, vt- The black solid lines 
show the strea klines resulting from the test particle s being released with the same angular velocity as the cluster. The epicyclic motion 
as described in lKiipper. MacLeod fc Heggij l l2008al ) is nicely reproduced. The red solid lines show the path of the test particles which 
have the same orbital velocity as the cluster. Just like the cluster they move on circular orbits about the galactic centre. The blue solid 
lines show the test particles with velocities in between these two cases. The length of the epicycles is shorter and the radial amplitude 
is smaller. The cluster mass was neglected in these calculations. Right panel: here the effect of the cluster mass on the trajectories of 
the test particles is shown. The cluster mass was chosen such that the tidal radius is 60 pc. The black solid lines are the same as in the 
left panel, whereas the red solid lines show the influence of the cluster mass: the sharp cusp is replaced by a smooth minimum, which, 
in addition, is further away from the cluster. Moreover, the test particles are, on average, further away from the cluster orbit in radial 
direction. 



bit about the galactic centre. Moreover, the cluster mass is 
neglected in this calculation. 

For eccentric cluster orbits such a general analytical so- 
lution does not exist. In addition, the trajectories of escaping 
stars depend on the orbital phase of the cluster in this case. 
Instead of solving the equations of motion analytically, we 
can simply illustrate the trajectories of stars for eccentric 
cluster orbits by using streaklines. This is a concept from 
fluid dynamics in which test particles are released into a 
fluid from a given point to visualise the flow. In engineering 
streaklines are often produced with smoke or dye. 

To produce the streaklines we integrate the cluster or- 
bit and, at given time intervals, release test particles. These 
particles have a certain velocity, v, and are released from two 
points which are given by the cluster centre plus a flxed pos- 
itive or negative radial offset. Ax, along the galactocentric 
radius vector. All released test particles are then integrated 
together with the cluster to the point in time of the snap- 
shot. The gravitational influence of the cluster on the test 
particles is taken into account in the cases where it is men- 
tioned explicitly. There is no interaction between the test 
particles. 

For this purpose, it is not reasonable to stay within the 
rotating, accelerated reference frame in which a cluster on a 
circular orbit is at rest, for which th e ana lytical solut i ons o f 
iKiipper. MacLeod fc Heggid l|2008al ) and lJust et alj (|2009l ) 
have been developed. Thus, we have to keep in mind that we 
have to convert the tangential offset velocities in the rotating 



coordinate system, vt, to the Cartesian, galactocentric rest 
frame velocities, v, when applying Eq. (6] 

First we illustrate that the technique works correctly by 
comparing it with results from the formalism described in 
Sec. l2.1l for circular orbits. Therefore, we produce three sets 
of streaklines in which the test particles are released with 
three different velocities, v, from the cluster on a circular 
orbit: 

(i) the test particles are released with the same angular 
velocity as the cluster, Q,, i.e. the test particles which are 
released from the larger galactocentric radius have a higher 
orbital velocity than the cluster, v = Vort + ilAx, whereas 
the test particles at the smaller galactocentric radius have a 
lower orbital velocity, v = Vorb — ^lAx. This corresponds to 
the case vt = in Eq. [G] 

(ii) the test particles have the same orbital velocity as 
the cluster, i.e. v = Vorb- This corresponds to vt = —QAx in 
Eq.El 

(iii) the test particles have an orbital velocity in between 
the two cases, i.e. vt = —O.bQAx corresponding in the galac- 
tocentric rest frame to u = Vorb + 0.5Q,Ax for the leading 
tail and v = Vorb — 0.5QAx for the trailing tail, respectively. 

These sets of streaklines we also produce for the eccentric 
orbits to see how the orbital acceleration and deceleration 
inffuence the trajectories of escaping stars. 

Finally, we produce additional sets of streaklines for 
which the mass of the cluster is taken into account. The 
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Figure 3. The same as in tiie left panel of Fig. [2] but for the cluster with an orbital eccentricity of 0.25. The four snapshots show the 
cluster and the corresponding streaklines in four different orbital phases (see also Tab.[TJ. Ax was chosen to be 60 pc just like in Fig. [2] 
The epicyclic movement gets distorted due to the eccentric cluster orbit. At apogalacticon (upper left panel) the epicycles are squeezed 
together, whereas at perigalacticon (upper right panel) the streaklines are stretched. In between these two orbital phases the shape of 
the streaklines differs if the cluster and its tails are being accelerated (lower left panel) or decelerated (lower right panel). Moreover, we 
get epicyclic movement also for the case when the test particles are released from the cluster with the same orbital velocity as the cluster 
(red solid lines). 



cluster mass has so far been neglected but will, of course, 
influence the trajectories of escaping stars. The cluster can 
attract stars and even re-capture them if they happen to lie 
beyon d the tidal radius (g iven by Eq.[2]) at some later orbital 
phase. iRead et al.l l|2006l ). for example, demonstrated that 
stars orbiting the cluster on prograde orbits with respect to 
the cluster motion about the galactic centre can escape the 
cluster potential more easily than stars on retrograde orbits. 
We therefore have to consider from which radius we release 
the test particles to produce the streaklines. 



If we chose to release the test particles from the 
apogalactic tidal radius, which is the largest tidal radius 
throughout the orbit since the cluster has the lowest angu- 
lar velocity at this point (see Tab.[l|, then we may expect 
that virtually all test particles will escape the gravitational 
attraction of the cluster mass. This, however, might not be 
the most realistic approach. If we, on the other hand, chose 
to release the test particles always from the actual tidal ra- 
dius, XL{t) (Eq.[2]), at the given time, t, then many test par- 
ticles may be re-captured. This will especially be a problem 



More on the structure of tidal tails 



Apogalacticon 





Between apo- 
and perigalacticon 



— I 1 1 1 

Same angular velocity, Q. 
Same orbital velocity, Vg^) 
Velocity in between 
Cluster orbit 



Perigalacticon 



Between peri- 
and apogalacticon 




-750 -500 



-250 250 
x[pc] 



500 750 



-750 -500 



-250 250 
x[pc] 



500 750 



Figure 4. The same as Fig. [3] but for an orbital eccentricity of 0.5. Again, the four snapshots show four different orbital phases (see also 
Tab. [ij. In this case Ax was also chosen to be 60 pc just like in Fig. [2] We see that the distortion effects get stronger with increasing 
eccentricity. The compression at apogalacticon is stronger than in Fig. |3] and so is the stretching at perigalacticon. Moreover, we can 
observe the effect of differential acceleration along the cluster-tail system. For example in the lower right panel the leading tail (to the 
right) is broader than the trailing tail. 



for the test particles released near perigalacticon since the 
tidal radius is smallest here and, in addition, changes rapidly 
in this orbital phase as the angular velocity is largest at this 
point. 

As we will demonstrate later, for this reason it is nec- 
essary to introduce a minimum radius from which test par- 
ticles can be released such that no test particles gets re- 
captured. This radius we call the "edge" radius, Xedge, and 
it serves as a lower limit to the actual radius from which test 
particles will be released, i.e. 

Ax{t) = max{XL{t), Xedge) ■ (15) 



We determine the value of this edge radius by incrementally 
increasing the lower limit from the value of the perigalactic 
tidal radius until no test particle is re-captured. 

In the following we will show that, for circular cluster 
orbits, this edge radius coincides with the constant tidal ra- 
dius, since any particle released from rest from a smaller 
cluster radius will not be able to escape from the cluster 
potential (for the given escape conditions |(i)|(iii)[ | . For ec- 
centric orbits this value is smaller than the apogalactic tidal 
radius but significantly larger than the perigalactic tidal ra- 
dius. Its value decreases with increasing eccentricity for fixed 
apocentric distance. 
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X [pc] X [pc] 

Figure 5. The same as Fig.[3]but here the cluster mass is taken into account. For the black solid lines the test particles were released 
from the cluster's apogalactic tidal radius. The cluster mass was chosen such that the apogalactic tidal radius is 60 pc. The lines are 
similar to the lines for which the cluster mass was neglected (black solid lines in Fig.fSj. When the test particles are released from the 
actual tidal radius of the cluster (red solid lines), the influence of the cluster mass gets more pronounced, as can be seen in the additional 
loops and bends. To produce these streaklines we had to introduce a minimum radius from which particles were released, Xf,(;gei since 
else some test particles would have been re-captured by the cluster at a later orbital phase. This "edge" radius is 49.0 pc, whereas its 
perigalactic tidal radius is 37.5 pc. 



2.4 Circular orbit 

In the left panel of Fig. [5] the resulting streaklines for the 
circular cluster orbit are shown. The dashed line gives the 
orbit of the cluster. Along this orbit test particles have been 
released in constant time intervals with a positive or neg- 
ative spatial offset, Aa;, from the cluster centre along the 
galactocentric radius vector and with one of the three veloc- 
ity offsets, Vt, given above. We chose Aa; to be 60 pc, which 
is a somewhat arbitrary choice. The cluster and the test par- 
ticles have each been independently integrated to the point 



in time of the snapshot. Hence, there is no gravitational in- 
fluence of the cluster on the test particles. 

From the left panel of Fig. it can be seen that th e an- 
alytical results o f iKiipper. MacLeod fc Heggid (|2008al ') and 
lJust et all l|2009l ) are nicely reproduced when the test parti- 
cles are released with the same angular velocity as the clus- 
ter (black solid line) . Moreover, if the orbital velocity of the 
test particles at the larger (smaller) galactocentric radius is 
lower (higher) than in the case of equal angular velocity, the 
length of the epicycle is shorter (blue solid line). 

In the case when the test particles have the same orbital 
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Figure 6. The same as Fig. [5] but for tiie cluster on an orbit witti eccentricity of 0.5. Again, the cluster mass was chosen such that the 
apogalactic tidal radius is 60 pc. The streaklines for which the test particles were released from the apogalactic tidal radius (black solid 
lines) are again similar to the streaklines for which the cluster mass was neglected (black solid lines in Fig.|4]l. When the test particles 
are released from the actual tidal radius the streaklines get much more complex (red solid lines). For these streaklines we had to assume 
an edge radius 32 pc, whereas the perigalactic tidal radius is 21.2 pc. 



velocity as the cluster then the epicycles should, in principle, 
be half as short as in the case of sharing the cluster's angular 
velocity, since 



yc ~ Stv [ Ax 



2Q 



= Stt- 



Ax 



(16) 



Instead, we see that the the test particles follow a circular 
orbit (red solid line). From Eq. [8] we see that, in the case 
when vt = —flAx, the radial extrema in the epicycles are at 



Vt 



x^ = 2Ax + ^ ± 



= Ax±Q. (17) 
Hence, the radial amplitude is zero, and the test particles 



are moving on a circular orbit about the galactic centre with 
a radial offset to the cluster of size Aa;. 

The theoretical framework (Eq. [3]-[8]) is in perfect agree- 
ment with our streaklines. Now we include the gravitational 
attraction of the cluster in our computations. The mass of 
the cluster in this computation is chosen such that the tidal 
radius, xl (Eq. [2}, coincides with the starting point of the 
test particles, Aa;, i.e. 60 pc in our example. 

As we can see in the right panel of Fig. (2] the cluster 
mass increases the length of the epicycles, yc, by about 25%. 
Moreover, the sharp cusp at yc is replaced by a smooth min- 
imum, and the radial distance of the test particles from the 
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cluster orbit, Ai, is larger on average. The radial amplitude 
appears to be damped by a factor of about 0.6 when the 
cluster mass is taken into account. 

All these effects can be explained when looking at the 
shape of the effective potential of cl uster and galaxy near 
the tidal radius (see e.g. figure 2 of I Just et aiTl2009h . Be- 



yond the tidal radius the effective potential is repulsive such 
that the te s t part icles get accelerated away from the cluster. 
ljust et al] (|2009|) found that, due to Jacobi-energy conser- 
vation (see e.g. Heggie fc Hut|[2003l ). a test particle starting 
at Aa; = xl with vt — will be at Aa; — \/Zxl if it ever 
comes to rest again in the future. In fact, the test particles 
in the minimum of our example are less than f 0% off this 
value as they still have a velocity of 0.5 km/s with respect 
to their initial conditions. 

Note that, since the distance of the Lagrange points as 
well as the length of the epicycles scales with the mass of 
the cluster, the effect of the cluster mass on the trajectories 
will always be the same. Stars escaping from a larger radius 
than the tidal radius will, of course, be less influenced by 
the cluster mass. Their trajectories will look more like the 
ones for which we neglected the cluster mass. 



2.5 Eccentric orbits 

In Fig. [3] the streaklines of the cluster in an eccentric orbit 
with eccentricity 0.25 are shown for four different orbital 
phases. Here, again, the cluster mass was neglected. Also, 
we chose Aa; to be 60 pc like in the previous example. We 
see that the streaklines are distorted in comparison to the 
circular orbit case (Fig.[2|. Moreover, the distortion depends 
on the orbital phase of the cluster. At apogalacticon the 
streaklines appear to be similar to the circular orbit case 
but squeezed together. The length of the epicycles is much 
shorter than in the circular orbit case and epicyclic loops 
further away from the cluster are tilted towards the cluster 
orbit. 

Between apo- and perigalacticon, when the cluster and 
the test particles are accelerated towards the galactic cen- 
tre, the tilt of the epicyclic loops gets stronger, whereas 
the length of the epicycles increases. At perigalacticon the 
epicyclic loops are almost completely aligned with the clus- 
ter orbit. The length of the epicycles is largest in this orbital 
phase. Between peri- and apogalacticon, when the cluster 
and the test particles are decelerated, the tilt gets weaker, 
the epicycles shorter and the loop structure gets more pro- 
nounced. 

Also apparent from Fig. [3] is the thickness of the tidal 
tails. Assuming that most stars follow such an epicyclic mo- 
tion, we can expect the tidal tails to be thickest at apogalac- 
ticon and thinnest at perigalacticon. Moreover, we see that 
we get epicyclic motion even when the escaping stars have 
the same orbital velocity as the cluster. This is due to the 
stars and the cluster being on slightly different eccentric or- 
bits about the galactic centre. The same is shown in Fig. [4] 
for the cluster with an orbital eccentricity of 0.5. All effects 
described above get more pronounced. Furthermore, with 
increasing orbital eccentricity it gets more difficult to sepa- 
rate the individual streaklines as the epicyclic loops overlap 
multiply. 

In these illustrations the offset radius. As, from which 
the test particles have been released were the same for all 



snapshots. When we take the cluster mass into account in 
our computations then we can fix Ax to the apogalactic tidal 
radius in order to be able to better compare the resulting 
streaklines to Fig. [3] & U Therefore we chose the cluster 
mass such that the apogalactic tidal radius is 60 pc. The 
streaklines are shown in Fig. [S] (black solid lines) for the 
cluster on an orbit with eccentricity of 0.25. As we can see, 
the streaklines look quite similar to the lines for which the 
cluster mass was neglected (black solid lines in Fig.jS]). Only 
some minor bends in the lines tell us that there is some 
influence of the cluster on the test particles. 

This changes when we release particles from the actual 
tidal radius (Eq. The influence of the cluster mass on 
the test particles gets much more obvious, as the lines show 
additional loops and bends (red solid lines). To produce the 
streaklines in Fig. [S] we had to introduce a minimum radius 
from which we release particles, the "edge" radius, Xedge, 
mentioned above (Sec. 12.31 Eq. II5p . For the orbit with an 
eccentricity of 0.25 this edge radius was found to be 49.0 
pc (about 80% of the apogalactic tidal radius), whereas its 
perigalactic tidal radius is 37.5 pc. 

Both sets of streaklines are shown in Fig.|S]for the clus- 
ter on an orbit with an eccentricity of 0.5. Here again we 
adjusted the cluster mass such that the apogalactic tidal 
radius is 60 pc. The effects described above get more pro- 
nounced. For the black solid lines the test particles were 
released from the apogalactic tidal radius, and for the red 
solid lines they were released from the actual tidal radius. 
The edge radius, which we had to introduce here as well, is 
32.0 pc (about 50% of the apogalactic tidal radius), whereas 
the perigalactic tidal radius is 2f .2 pc. 



3 COMPARISON WITH iV-BODY DATA 

The three A''-body computations which we use here for 
compar ison with streaklines are taken from iKiipper et al.l 
lioioi). The models were set-up using the publicly avail- 
able code McLustefQ l|Kiipper et al.ll2011bl ) . They initially 
consist of 65536 stars d rawn from the canonical initial mass 
function (|Kroupall200"ll . eq. 2) with masses between 0.1 M© 
and 1.2 M©, resulting in a cluster mass of about 20000 M0. 
The cluster stars follow a Plummer density profile with a 
half-mass radius of 8 pc. 

The clusters are initially located at their apoc- 
entre of 8.5 kpc with in an Milky Way-like potential 
l| Allen fc SantillanllT99ll : see Sec. 12. 2p . We chose three or- 
bital eccentricities for the clusters of e = {0.00,0.25,0.50}, 
which start with an orbital velocity of Vorb ~ {220, 165, 110} 
kms"'^. From Eq.[T]we can infer that the tidal radius of these 
clusters is of the order of 40 pc. Hence, with a ratio of half- 
mass radius to tidal radius of about 0.2 the models belong to 
th e group of extended, tida lly filling clusters as categorised 
bv lBaumgardt et al.l (I2OI0I ). 

The computations were performed with Nbody4 
( AarsethI 20031) on the Grape sp ecial-purpose computers 



( Fukushige. Makino fc KawailbOOSl ) at AlfA Bonn for 4 Gyr. 
Here we use snapshots in the middle of these computations 



^ www . astro . uni-bonn . de/'akuepper/mcluster/mcluster . html 
or www. astro .uni-bonn. de/~webaiub/gerinan/downloads .php 
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Figure 7. Comparison of streaklines to results from A'^-body com- 
putations of a cluster on a circular orbit. The cluster was inte- 
grated for 2 Gyr and has lost about 3600 Mq of stars which arc 
mainly in the tidal tails now (black dots; stars within a cluster 
radius of 40 pc were omitted for clarity). The blue solid lines 
show the streaklines for which the cluster mass was neglected. To 
match the length of the epicyclic overdensities, the test particles 
were released with a tangential velocity of vt = 0.55 kms~^ from 
the cluster's tidal radius, i.e. Ax = xj^ = 37 pc. For the red solid 
lines the cluster mass was taken into account. The test particles 
were released from the cluster's tidal radius with zero velocity. 
As we can see, the shape of the tails is much better reproduced 
by the streaklines which take the cluster mass into account. The 
scatter of stars about this streakline may well originate from the 
scatter in escape conditions. 



at different orbital phases as shown in Fig.[l]and summarised 
in Tab. [T] 



3.1 Circular orbit 

The first snapshot (Fig. [7]) is from the cluster on a circular 
orbit after 2 Gyr of dynamical evolution. The star cluster is 
shown in projection onto the orbital plane similar to Fig. [J] 
It has lost about 3600 Mq , of which most stars escaped into 
the tidal tails (black dots in Fig. [T] stars within a cluster 
radius of 40 pc are not shown for clarity). The predicted 
over- and underdensity pattern is clearly visible in both the 
leading and the trailing tail. Moreover, we see that the tidal 
tails have a rolling shape, resulting from the epicyclic motion 
of the tail stars. 

The tidal radius, xl (Eq. [2]), of the cluster at this mo- 
ment is 37 pc, so from Eq. |4] we expect the length of the 
epicycles, yc, to be about 350 pc when the stars escape 
with zero tangential velocit y, Vt- We find, howeve r, that they 
are closer to 400-500 pc. In lKiipper et al.l (j2010ah we argued 
that this is due to the cluster stars escaping with small excess 
velocities of about 0.27-0.55 kms~^ in the direction of the 
tails. This is a reasonable assumption, since such excess ve- 
locities have been measured for evaporating s t ars in A '^-body 
simulations. iKiipper. Kroupa fc Baumgardtl l|2008bl ) found 
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Figure 8. The same as Fig. [7] but in velocity space (velocities 
arc given with respect to the cluster velocity). Black dots show 
more or less the same stars from the A^-body simulation as shown 
in Fig. [T] Here stars within a cluster radius of 40 pc were again 
omitted for clarity. As in Fig. [7] the blue solid lines show the 
streaklines for which the cluster mass was neglected, whereas for 
the red solid lines the cluster mass was taken into account. Also 
in velocity space the shape of the tails is much better reproduced 
by the streaklines which take the cluster mass into account. 



the velocities of evaporating stars to be of a log-normal dis- 
tribution with a peak at roughly the mean stellar velocity 
within the cluster. Also, it is clear from Fig. [2] that, in our 
simple model, the distance between overdensities depends 
on the velocity with which escapers are released. Neverthe- 
less the streaklines shown in Fig. [7] were generated from a 
single value of the initial velocity, which we have found to 
produce a satisfactory fit to the structure of the tidal tails. 

The streaklines of these assumed escape conditions are 
shown in blue in Fig. [T] As we can see, the epicyclic pat- 
tern is well reproduced, that is the epicycles have the cor- 
rect length. Nevertheless, these streaklines do not seem to 
perfectly describe the motion of escaping stars. It seems like 
most escapers rather follow trajectories which are less curved 
and which are closer to the cluster orbit on average. How- 
ever, we neglected the cluster mass for these streaklines. 

For the red solid lines in Fig. [7] the cluster mass was 
taken into account. The test particles were released with zero 
velocity from the tidal radius in this case. The streaklines 
seem to perfectly match the mean path of the escaping stars, 
therefore they refiect the dominant mode of escape. The 
scatter about this mean path may well be due to the scatter 
in escape conditions mentioned above. 

In Fig.[8]the comparison between the A''-body data and 
our simple streakline model is shown as in Fig. [7] but now 
in velocity space. The epicyclic motion of tail stars can also 
be seen as a periodic velocity variation. Like in Fig. [7] we 
can see that the A'^-body data is reproduced better by the 
streaklines for which the cluster mass is taken into account. 
Our simple model therefore seems to represent the average 
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escape conditions. As can be seen in Figs. [7| & [S] the de- 
viations from these average escape conditions are relatively 
small. 



3.2 Eccentric orbits 

In Fig.[5]the four consecutive snapshots of the cluster on an 
orbit with an eccentricity of 0.25 are shown. The cluster has 
lost about 5000 of stars which are mainly in the tidal 
tails (black dots). Here, again, the epicyclic movement of 
the tail stars is visible. The tails show periodic overdensities 
with a distance from th e cluster which varies with the orbital 
period (as observed bv lKiipper et aLlboiOal ). At apogalacti- 
con the overdensities are close to the cluster, whereas they 
are furthest away at perigalacticon. 

Also the shape of the tails depends on the orbital phase. 
At apogalacticon they appear similar to the circular orbit 
case with the rolling shape. However, the tails are denser 
and broader as they are being compressed due to differential 
acceleration across the tails. At perigalacticon the tails are 
thinnest and they appear almost like straight lines. The two 
snapshots in between appear quite different even though the 
cluster is in a similar orbital phase, only being accelerated in 
the one snapshot and decelerated in the other. The distance 
of the overdensities is comparable, whereas the shape of the 
tails is very distinct. In between peri- and apogalacticon 
the tails show prominent streaks like ha ve often been 
observed in numerical simulations before (iDehnen et al] 



2004 'Capuzzo Dolcetta. Pi Matteo fc Miocchil l2005l : 
Montuori ct al. 2003)5 but which have not yet been assigned 



to epicyclic motion of tail stars. These features may also 
be the true nature of the multiple tidal tails observed for 
Galac tic satellites like NGC 288, Willman 1 and NGC 
2298 ('Loon. Mcvl an fc Combed l2000l : I Willman etal] l2006l : 
[Balbinot ct al. 20l"ll '). 

Also shown are streaklines for which the test particles 
have been released from the actual tidal radius (Eq.[2|. How- 
ever, since the tidal radius at perigalacticon is very small 
(24.8 pc) and the passage through perigalacticon is quite 
fast, there are many test particles which get re-captured 
from the cluster as it moves to apogalacticon. Therefore, we 
had to introduce a minimum radius, x^dge, from which test 
particles are released (see Sec. 12.31 Eg. I15|) . This "edge" ra- 
dius was found to be 32.5 pc by incrementing its size from 
zero until no test particles was re-captured. Its size is about 
80% of the cluster's apogalactic tidal radius (39.5 pc). Re- 
markably, with this minimal choice of edge radius the streak- 
lines reproduce the TV-body data very well, even though their 
shape is very sensitive to this value. 

In Fig. [To] the four consecutive snapshots of the cluster 
on an orbit with e = 0.5 are shown. This cluster has already 
lost 10500 M© of stars which are now mainly in the tidal 
tails (black dots). The orbital compression of the tails gets 
stronger at apogalacticon. The overdensities are very close 
to each other and the tails are very broad. Moreover, they 
show prominent streaky features. At perigalacticon the tails 
become even thinner than in the e = 0.25 case. From the 
snapshot between peri- and apogalacticon we can see that 
the cluster-tail system is experiencing severe differential ac- 
celeration, such that the leading tail (to the right) is already 
being broadened, whereas the trailing tail is still quite thin. 

The red solid lines in Fig. llOlhave been produced by re- 



leasing test particles from the actual tidal radius and assum- 
ing a minimum tidal radius (edge radius) of 18.0 pc fEg. fTS)) . 
The perigalactic tidal radius of this cluster is 12.9 pc and the 
apogalactic tidal radius is 36.4 pc. Hence, the edge radius 
is about 50% of the apogalactic tidal radius. Again, as for 
the cluster with an orbital eccentricity of 0.25, this choice of 
escape conditions reproduces the A'^-body data remarkably 
well. 

Finally, from Fig. [9] & [10] it is obvious that all substruc- 
ture within the tidal tails can be ascribed to epicyclic motion 
of tail stars. This rules out tidal shocks as origin of substruc- 
ture in these systems. Moreover, it demonstrates how com- 
plex the motion of stars within tidal tails of star clusters 
like Pal 5 can be, and thus how unlikely the formation of 
overdensities in such tidal tai l s through Jeans instabiliti es is 
(iQuillen fc Comparettall2010l : [Schneider fc Moordl201ll ). 



4 SUMMARY AND CONCLUSIONS 

We have investigated the motion of stars escaping from 
star clusters using streaklines (for details on streaklines see 
Sec. I2.3|l , focussing on three star cluster models with orbital 
eccentricities of 0.0, 0.25 and 0.5 (further description of the 
models in Sec. 12.21) . We demonstrated how escaping stars 
move on epicycles within tidal tails (Sec. l2.4l fc l2.5p , and that 
this epicyclic motion is the only origin of the substructure 
observed in Ai'-body computations of dissolving star clusters 
(Sec.[3|^ 

In [Kiipper. MacLeod fc Heggid l|2008al ) it was shown 
analytically and numerically for circular cluster orbits, how 
the length of the epicycles changes with the size of the 
radial offset, /S.x, escaping stars have with respect to the 
cluster orbit at the moment of escape. Here we visualised, 
with the help of streaklines, how the shape of the epicycles 
changes when we vary the tangential velocity, vt, of the es- 
caping stars. In fact, the shape of the epicycles and thus 
the tidal tails is very sensitive to this velocity. Furthermore, 
we demonstrated how this motion gets distorted when the 
cluster orbit is eccentric, and what (important) influence the 
cluster mass has. 

When a cluster on an eccentric orbit is close 
to apogalacticon then the epicyclic motion of stars 
within its tails may lead to streaks appearing 
like multiple tails. Such features hav e been ob- 
served i n numerical simulations before ([Pehnen et al] 



,200- 



Capuzzo Dolcetta. Pi Matteo fc Miocchil [2005[ : 



Montuori et al 



2007[ ) 



but have not yet been interpreted 
as result of epicyclic motion. The streaks may also be 
the true nature of the multiple tidal tails observed for 
Galac tic satellites like NGC 288, Willman 1 and NGC 
2298 ([Leon. Mevlaii fc Combed [2000l : [Willman etall [2006[ : 
[Balbinot et al.|[201l[ '). 

Finally, we compared sets of streaklines to results from 
three A'^-body models following the same orbits as stated 
above. We found that the tidal tails of these computations 
can be accurately reproduced by a quite simple model, that 
is when we assume that the stars evaporate from the actual 
tidal radius, XL{t), the star cluster has at the moment when 
the star escapes. When the velocity of the escapers is chosen 
such that they have the same angular velocity as the cluster. 
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Figure 9. The same as Fig. [7] but for the cluster with an orbital eccentricity of 0.25. This cluster has lost already 5000 of stars 
(black dots; stars within a cluster radius of 40 pc were omitted for clarity). As in Fig. [3]&[5] the four snapshots show the cluster in 
four different orbital phases. The shape and density of the tidal tails changes significantly during one orbit. At apogalacticon they look 
similar to the circular orbit case but squeezed together. At perigalacticon the tails are stretched and the overdensities are much further 
away from the cluster than at apogalacticon. The orbital phases in between look quite distinct: between peri- and apogalacticon the tails 
show a streaky structure which has often been observed in numerical simulations. Also shown are streaklines for which the test particles 
have been released from the actual tidal radius (Eg. HSU with zero velocity (red solid lines). The edge radius which was used here is 32.5 
pc, in comparison to the 24.8 pc perigalactic tidal radius and the 39.5 pc apogalactic tidal radius. This set-up can reproduce all four 
snapshots without further ado, and especially reproduces the streaky features between peri- and apogalacticon. Such features have often 
been interpreted as results of tidal shocks in the past. 



the shape of the tidal tails is reproduced best. This is what 
we find to be the dominant mode of escape. 

However, we have to assume that there is a lower 
limit, Xedge, to the tidal radius. That is, stars escape from 
Ax{t) = max{xLit),Xedge)- This "edge" radius was found 
to be identical to the tidal radius of the cluster, if it moves 
on a circular orbit, but was found to decrease with increas- 
ing orbital eccentricity at fixed apocentric distance. For all 



orbits, the edge radius turned out to be significantly larger 
than the corresponding perigalactic tidal radius. It can be 
interpreted as the minimum radius from which a star can es- 
cape from the cluster during one cluster orbit without being 
re-captured by the cluster in a later orbital phase. 

These findings have only been established, however, for 
systems where the escape of stars is dominated by collisional 
heating, on galactic orbits with eccentricities up to 0.5. For 



14 A.H.W. Kiipper, R.R. Lane, B.C. Heggie 



750 



500 



250 



o 

a 



-250 
-500 
-750 

750 
500 
250 



o 

S 



-250 
-500 
-750 





Apogalacticon 




Between apo- 
and perigalacticon 



1 1 1 1 — 

N-body + 
Ax = max(XL(t), Xedge) 



Perigalacticon 



H h 



Between peri- 
and apogalacticon 




-750 -500 



-250 250 
x[pc] 



500 750 



-750 -500 



-250 250 
x[pc] 



500 750 



Figure 10. The same as Fig. [7] & |9] but for an orbital eccentricity of 0.5. Due to the high eccentricity, this cluster has already lost 
10500 Mq of stars (black dots; stars within a cluster radius of 40 pc were omitted for clarity). Here, like in Fig. |4l the four snapshots 
show four different orbital phases. We see that the compression and stretching of the tails gets more pronounced with increasing orbital 
eccentricity. Furthermore, at apogalacticon the compressed epicyclic loops appear like streaky features. These structures can be well 
reproduced by streaklines, when we release the the test particles from the actual tidal radius (red solid lines). The edge radius wc used 
here is 18.0 pc, in comparison to the 12.9 pc perigalactic tidal radius and the 36.4 pc apogalactic tidal radius. 



much larger eccentricities, or any other system where mass 
loss may be dominated by tidal heating, modelling the struc- 
ture of tidal tails may requi re refinement, as sugges ted by 
the collisionless modelling in lPenarrubia et al.l (|2009l ). Even 
in the case of a cluster on a circular galactic or bit, the no- 
tion of a single escape radius is over-simplistic (|Read et al.l 
but even these authors find that, in the long run, the 
"tidal radii" of stars on prograde, circular and retrograde 
orbits tend to a single value. 

iKiipper et al] (|2010br ) investigated in detail the same 
three A'^-body models that we discuss in the present pa- 



per. They found that all three models show a radial den- 
sity pro file wh ich barely changes with time. By fitting King 
profiles jKind [196211 ) and a new, so-called, KKBH profile, 
they showed that the fitted tidal radius does not change 
with the orbital period, even though Eq. [2] would suggest 
that the tidal radius should be large at apogalacticon and 
small at perigalacticon. For this reason they suggested to 
rather name such a fitted tidal radius the "edge" radius. 
Furthermore, they found that the fitted edge radius of the 
model on a circular orbit is close to its true tidal radius 
of 37 pc, whereas the model with an orbital eccentricity of 
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0.25 showed an edge radius of about 30 pc, and the one with 
e = 0.5 had a constant edge radius of about 20 pc. These val- 
ues he remarkably close to the values we found by adjusting 
streaklines to the tidal tails of these A'^-body models. 

These results suggest that the edge radius is fundamen- 
tally connected to the orbit of a star cluster and its mass. 
Stars preferentially evaporate from the actual tidal radius 
of the cluster whose lower limit is given by the edge radius. 
That is, stars within the edge radius cannot escape, whereas 
stars outside the edge radius get "eaten away" by the tide. 
This edge radius can be recovered by fitting a King or a 
KKBIf profile to the radial density profile. It must be men- 
tioned, though, that these findings will mo st likely only hold 
for e xtended, tidally filling clusters (see iBaumgardt et all 
with a ratio of half-mass radius to tidal radius larger 
than approximately O.I, since more concentrated clusters 
may not have properly adapted to the mean tidal field yet. 

Furthermore, applying our results to dwarf spheroidal 
galaxies may be complicated by the fact that in such systems 
the crossing time of stars orbiting at large radii may easily 
exceed the orbital time of the dwarf galaxy about its host. 
Instead, Pcfiarru bia et al.l ((2009, ) find that in dwarf galax- 
ies which experienced a recent tidal shock an edge devel- 
ops which moves radially outward with time. Such a dwarf 
galaxy will lose most of its mass due to tidal stripping rather 
than through dynamical evaporation. The scatter in escape 
conditions of escaping stars would be much higher in this 
case so that our simple model would not work for the ma- 
jority of escapers. Thus, no pronounced epicyclic signature 
would be expe cted for the tidal tails of such objects. 

Originally, lKing|l|l962l ) suggested that the observed cut- 
off radius, i.e. edge radius, of a star cluster on an eccentric 
cluster orbit is related to its respective perigalactic tidal 
radius 

where f2peri is the cluster's angular velocity at perigalacti- 
con. Hence, the cut-off radius should be even significantly 
smaller than the perigalactic tidal radius for high eccentric - 
ities . Similar results have been found by I Read" et all (|2006l ') 
and iKenne d3 (l20Ilh . both predicting the cut-off radius to 
be of the order of the perigalactic tidal radius. In contrast 
to this, we found the edge radius of clusters on eccentric 
orbits to be significantly larger than the perigalactic radius. 
The reason for this difference is that stars that become un- 
bound at perigalacticon at smaller radii get re-captured by 
the cluster at a later orbital phase. 

Finally we conclude that, if sufficient information on a 
specific globular cluster of the Milky Way is available, this 
knowledge on the formation of tidal debris can be used to 
predict the shape of its tidal tails and the location of its 
epicyclic overdensities. In a future investigation we are go- 
ing to demonstrate this with the globular cluster 47 Tucanae 
(Lane, Kiipper & Heggie, in prep.). Moreover, we would like 
to point out that epicyclic motion of stars within tidal tails 
of star clusters leads to over- and underdensities indepen- 
dent from the cluster orbit about the host galaxy. Therefore 
substructure in tidal tails is not necessarily due to substruc- 
ture in the galactic potential like spiral ar ms, giant molecu- 
lar clouds or d ark-matter subhaloes fe.g. Ilbata et al.l[2002l : 
ICarlberdlioogl ). 
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